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' Abstract 



The inclusive production cross sections for W + , W~ and Z°-bosons form important bench- 
marks for the physics at hadron colliders. We perform a detailed comparison of the pre- 
dictions for these standard candles based on recent next-to-next-to-leading order (NNLO) 
parton parameterizations and new analyses including the combined HERA data, compare 
to all available experimental results, and discuss the predictions for present and upcoming 
RHIC, SPS, Tevatron and LHC energies. The rates for gauge boson production at the LHC 
can be rather confidently predicted with an accuracy of better than about 10% at NNLO. 
We also present detailed NNLO predictions for the Higgs boson production cross sections 
for Tevatron and LHC energies (1.96, 7, 8, 14 TeV), and propose a possible method to mon- 
itor the gluon distribution experimentally in the kinematic region close to the mass range 
expected for the Higgs boson. The production cross sections of the Higgs boson at the LHC 
are presently predicted with an accuracy of about 10-17%. The inclusion of the NNLO 
contributions is mandatory for achieving such accuracies since the total uncertainties are 
substantially larger at NLO. 
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1 Introduction. The inclusive production cross sections of the weak gauge bosons W ± ,Z° 
and of the Higgs particle H° form important reference points for the physics at hadron colliders. 
Within QCD the corresponding cross sections have been calculated to NNLO [1,2]. 2 This 
level of accuracy is necessary to control the renormalization and factorization scale uncertainties 
of the parton distribution functions (pdfs), which are still of significant size at next-to-leading 
order (NLO). The rapidly growing luminosity at the Large Hadron Collider, LHC, allows a precise 
measurement of these quantities [6,7], which reaches the accuracy of the NNLO predictions, based 
on the world deep-inelastic, Drell-Yan and di-muon data, and part of the Tevatron data [8-13]. 
These analyses also require a correct treatment of the heavy flavor contributions, cf. [14], which 
yield sizeable effects. 

In this note we provide NNLO predictions for the inclusive weak boson and Standard Model 
Higgs-boson production cross sections, and compare to the available experimental data [6,7,15- 
20]. For present and upcoming experimental analyses we also provide detailed reference tables, 
and comment on the well-known NLO predictions, cf. [21]. These are, however, far less precise 
due to their inherent large renormalization and factorization scale uncertainties [22,23]. 

The effect of the new combined HERA data on the inclusive NNLO W ± , Z° and H° boson 
cross sections [24] is shown by comparing the results for the recent ABM10 [8] and ABKM09 [9] 
distributions. It leads to significant shifts of the total rate of about la in the pdf-error for a wide 
range of collider energies and processes. We consider all collider energies having been probed 
so far, compare to all measurements, and give predictions for the high energy options at the 
LHC. In this way a wide kinematic region in Bjorken x is probed for the corresponding parton 
luminosities which may help to delineate remaining differences in the current NNLO pdfs and 
to devise a way of potential further improvements. 

2 W ± and Z° Boson Production. The inclusive W ± and Z° boson production cross sections 
in pp and pp scattering are known to 2nd order (NNLO) in the strong coupling constant [1], 
supplemented by the 1st order (NLO) electroweak corrections, cf. [3]. In the following comparison 
we will concentrate on the QCD corrections only. The electroweak parameters are calculated 
choosing the scheme based on (G F ,M W ,M Z ) [25] with G F = 1.16637 x 10~ 5 GeV~ 2 , M w = 
80.399 ± 0.023 GeV, M z = 91.1876 ± 0.0021 GeV, and the weak mixing angle as dependent 
quantity, with 

M 2 

§1 = 1 - — ^ = 0.2307 ±0.0005 , (1) 

z pM\ v ' 

and p = 1.01047 ± 0.00015. Furthermore, the width of the W ± and Z° bosons are r(VF ± ) = 
2.085±0.042 GeV, T(Z°) = 2.4952±0.0023 GeV, and sin 2 6 C = 0.051, with 6 C the Cabibbo angle. 
We compute the inclusive production cross sections at various collider energies for the recent 
NNLO parton distributions, ABM10 [8], ABKM09 [9], JR [10], MSTW08 [11], and HERAPDF 
[12]. In the latter case we refer to the fit where a value of a s (M z ) = 0.1145 has been assumed 
at NNLO, which resulted in the lowest value for Xmm among other choices. 3 

In Table 1 we summarize the cross sections for the different pp collision energies at the SPS 
and the Tevatron at NNLO for the distributions [8-12] . For one of the parameterizations, JR, 
we compare also to the NLO QCD corrections. The relative differences between the NLO and 
NNLO corrections are of comparable size for all other pdf sets. We also list the corresponding 
NNLO values of a s (M|) and their errors as determined, or partly being assumed, in the various 

2 Electro-weak corrections were calculated in [3-5]. 

3 For the time being the NNLO HERAPDF parameterization yields only central values, unlike the NLO 
distributions, where also errors are provided, cf. [12]. 
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analyses, which in most cases turn out to be similar. We note that the corresponding a s values 
at NLO, despite differences in the central values, necessarily all agree within the rather large 
theory errors due to factorization and renormalization scale uncertainties of at least ±0.005, 
cf. [22,23]. This also applies to other NLO analyses [26,27]. 



(TeV) 




0.546 


0.630 


1.8 


1.96 


ABM10 [8] 


W ± 


5.632 ± 0.092 


7.045 ± 0.111 


24.441 ± 0.235 


26.740 ± 0.259 


a s = 0.1147 ±0.0012 


z° 


1.761 ± 0.022 


2.187 ± 0.028 


7.181 ± 0.068 


7.846 ± 0.075 


ABKM09 [9] 


w ± 


5.804 ± 0.075 


7.222 ± 0.091 


23.88 ± 0.243 


26.09 ± 0.265 


a s = 0.1135 ±0.0014 


z° 


1.806 ± 0.020 


2.234 ± 0.024 


7.056 ± 0.068 


7.691 ± 0.075 


JR [10] 


w ± 


5.983 ± 0.148 


7.346 ± 0.159 


23.069 ± 0.238 


25.157 ± 0.251 


a s = 0.1124 ±0.0020 




(5.358 ± 0.152) 


(6.637 ± 0.167) 


(22.121 ± 0.274) 


(24.181 ± 0.296) 




z° 


1.837 ± 0.029 


2.268 ± 0.034 


6.975 ± 0.071 


7.586 ± 0.076 






(1.648 ± 0.028) 


(2.047 ± 0.033) 


(6.667 ± 0.080) 


(7.272 ± 0.087) 


MSTW08 [11] 


w ± 


5.469 ± 0.151 


6.802 ± 0.176 


23.14 ± 0.394 


25.35 ± 0.422 


a s = 0.1171 ±0.014 


z° 


1.654 ± 0.047 


2.056 ± 0.056 


6.773 ± 0.126 


7.406 ± 0.134 


HERAPDF [12] 


w ± 


6.121 


7.519 


24.51 


26.80 


a s = 0.1145 


z° 


1.853 


2.296 


7.319 


7.978 



Table 1: NNLO predictions for the production cross sections a(pp — > V ± X) [nb], with V = W^, Z°. 
The abbreviation refers to the sum W + ± W~. Notice that for pp collisions the W + and W~ 
cross sections are equal. The errors refer to the ±lcx pdf uncertainties. The NNLO values of a s refer 
to a s = a s (Mz). To allow for a comparison with the corrections up to NLO the corresponding cross 
sections for the JR distributions are also listed as an example in parentheses [28]. 

The NNLO corrections furthermore enhance the cross sections at all center-of-mass (cms) 
energies both for W ± and Z° boson production. The importance of full NNLO analyses becomes 
evident when comparing the NLO and NNLO predictions in Tables l and 2. At the Tevatron 
(1.96 TeV), for example, Aa w± = a^ LO - a^o = ±0.930, ±0.976 and ±0.825 nb for ABM10, 
JR and MSTW08, respectively. Similarly for LHC (7 TeV) one obtains Aa w± = ±1.05, ±2.36 
and ±2.46 nb for ABM10, JR and MSTW08, respectively. This corresponds at the Tevatron to 
a difference of more than 3a w.r.t. the pdf-errors, and to more than la at LHC. Furthermore, at 
NNLO the parton distributions are much more stable against renormalization and factorization 
scale uncertainties than at NLO : the scale variations of the full NNLO predictions in Table 1 
amount to less than 0.5 % (i.e., are about half as large as the stated l<r p df uncertainties) which 
is about four times smaller than at NLO [10,28]; at LHC energies the scale uncertainties of the 
NNLO predictions in Table 2 amount to less than 2 % of the total predicted rates which is about 
half as large as the stated lcr p df uncertainties and the scale uncertainties at NLO [10,28]. 

Thus far, the combined HI and ZEUS data [24] have been taken into account in the ABM10 [8] 
and HERAPDF [12] analyses only, and their relative effect can be seen by comparing the numbers 
for the ABM10 and ABKM09 distributions at NNLO. While for the lower collider energies 2- 
3.5% lower cross sections are obtained, about 2.5 % larger cross sections result for the Tevatron 
energies of 1.8 and 1.96 TeV. With respect to the current pdf-errors this effect amounts to 2a. 
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y/s (TcV) 




0.5 


7 


10 


14 


ABM10 [8] 


W+ 


1.236 ± 0.057 


59.86 ± 0.838 


85.58 ± 1.267 


118.4 ± 1.891 




w~ 


0.363 ± 0.092 


40.28 ± 0.535 


60.28 ± 0.852 


86.58 ± 1.331 




w ± 


1.600 ± 0.070 


100.1 ± 1.315 


145.9 ± 2.065 


205.0 ± 3.186 




z° 


0.305 ± 0.015 


29.01 ± 0.391 


42.77 ± 0.633 


60.69 ± 0.963 


ABKM09 [9] 


W+ 


1.160 ± 0.046 


58.86 ± 0.903 


85.14 ± 1.427 


119.4 ± 2.072 




w~ 


0.348 ± 0.014 


39.43 ± 0.614 


59.56 ± 0.993 


86.53 ± 1.525 




w ± 


1.509 ± 0.058 


98.27 ± 1.527 


144.7 ± 2.436 


205.9 ± 3.658 




z° 


0.287 ± 0.012 


28.42 ± 0.457 


42.28 ± 0.743 


60.70 ± 0.115 


JR [10] 


W+ 


1.138 ± 0.061 


54.57 ± 1.10 


78.43 ± 1.98 


109.31 ± 3.13 






(1.245 ± 0.065) 


(52.96 ± 0.99) 


(76.60 ± 1.74) 


(107.58 ± 2.95) 




w~ 


0.387 ± 0.028 


37.15 ± 0.79 


55.54 ± 1.44 


80.02 ± 2.31 






(0.427 ± 0.030) 


(36.39 ± 0.72) 


(54.67 ± 1.26) 


(79.16 ± 2.12) 




w ± 


1.525 ± 0.052 


91.72 ± 1.82 


133.99 ± 3.35 


189.29 ± 5.41 






(1.672 ± 0.053) 


(89.36 ± 1.57) 


(131.23 ± 2.87) 


(186.74 ± 4.95) 






0.300 ± 0.011 


27.24 ± 0.50 


40.39 ± 0.95 


57.85 ± 1.56 






(0.336 ± 0.012) 


(26.57 ± 0.43) 


(39.57 ± 0.81) 


(57.00 ± 1.42) 


MSTW08 [11,13] 


W+ 


1.221 ± 0.0421 


56.80 ± 0.971 


81.83 ± 1.405 


114.0 ± 1.945 




w~ 


0.416 ± 0.017 


39.63 ± 0.678 


59.45 ± 1.008 


85.63 ± 1.484 




w ± 


1.637 ± 0.052 


96.41 ± 1.607 


141.3 ± 2.372 


199.6 ± 3.379 




Z° 


0.319 ± 0.011 


27.89 ± 0.481 


41.34 ± 0.705 


58.99 ± 1.012 


HERAPDF [12] 


W+ 


1.219 


59.37 


85.37 


119.0 




w~ 


0.414 


40.82 


61.06 


87.94 




w ± 


1.633 


100.2 


146.4 


206.9 




z° 


0.322 


29.08 


42.95 


61.22 



Table 2: NNLO predictions for the production cross sections a(pp — > V + X) [nb], with V = W ± , Z°. 
The abbreviation denotes the sum W + + W~ . The errors refer to the ±lcr pdf uncertainties. To 
allow for a comparison with the corrections up to NLO we also listed the corresponding cross sections 
for the JR distributions as an example in parentheses [28] . 

For illustration we compare in Figure 1(a) the different predictions [8,10-12] with each other 
and with the measured cross sections [15-20] at fixed energies in the range = 0.5 - 1.96 
TeV, covering also the pp cross sections at RHIC, cf. Table 2. A detailed summary of the 
various measurements is given in the Appendix. As can also be seen in Table 1 the results 
for Z° production predicted by the different parameterizations agree somewhat better than in 
the case of IT ± production. The experimental errors at lower energies v 7 ^ < 0.63 TeV are 
rather large and the current NNLO predictions agree with experiment. For v^S* = 0.63 TeV 
all NNLO analyses predict cross sections which are at the upper end of the measurements. At 
Tevatron energies the ABM10 and HERAPDF distributions yield about 5-7% larger W ± cross 
sections than those of JR and MSTW08. Moving to higher energies the JR distributions result 
in lower cross sections as compared to the ABM 10 and HERAPDF distributions, which yield 
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larger values at lower energies as well. Yet all predictions are in agreement with the current 
experimental results at the level of 1.5a. 



S 
\ 

£ 25 



10 - 



NNLO 

ABM10 

JR 

HERAPDF 
MSTW08 



tV*+MT , f 



o IM1 
□ UA2 



PHENIXW* 
PHENIXW 
CDF 
DO 



" I 1 "V 



:::: 



..i V 



to 



K 
S 



150 - 



125 - 



100 



75 



50 



25 - 



NNLO 

ABM10 

JR 

HERAPDF 
MSTW08 

ATLAS 
CMS 



W' + W 



W* 



W 

Z° 



* t 

t 


• * 


■1 . .. 




I 


■ 


■ 1? ■ 




,,iT 

•' 4 V 





(b) 



0.500 TeV 0.540 TeV 0.630 TeV 1.8 TeV 1.98 TeV 



7 TeV 10 TeV 14 TeV 



Figure 1: Comparison of different NNLO predictions for the inclusive W + , W , W^, and Z° boson 
production cross sections in pp annihilation and pp scattering (\fS = 0.5 TeV) based on the pdfs of 
recent NNLO analyses [8-13] and the corresponding experimental data [6,7,15-20]. Left panel (a): the 
lower energy region corresponds to pp collisions, except at 0.5 TeV, which refers to pp scattering. For 
the latter case the predictions refer to (from above) W + + W~ , W + , W~ and the ones for Z° are given 
to the right of the ones for W~ . Right panel (b): LHC energies {pp collisions); the inner error bars 
refer to (Cg tat + c^yst) 1 ^ and the total error is obtained by adding the luminosity error in quadrature. 

In Table 2 we list the NNLO predictions for the pp production cross sections for W + , W~ 
bosons, their sum, as well as those for Z° bosons at RHIC and for different present or planned 
collider energies at the LHC. As before we supplement the NLO prediction for one set of distribu- 
tions (JR) to allow for comparisons as an example. For the LHC energies the NNLO corrections 
lead to an enhancement of the cross sections w.r.t. NLO, while at RHIC energies the NNLO 
cross sections become smaller. At LHC energies the pattern of relative differences of the various 
predictions [8-13] both for W ± and Z° boson production pertains over the whole energy range, 
see also Figure 1(b). The cross sections grow empirically nearly linearly with ^/S. The impact 
of the combined HERA data is less than at lower collider energies and leads in most cases only 
to a very slight enhancement of the cross sections, as a comparison of the values in Table 2 re- 
sulting from the ABM10 and ABKM09 distributions shows. The largest scattering cross sections 
are obtained for the ABM10 and HERAPDF distributions, followed by slightly lower values for 
MSTW08, and even somewhat lower cross sections are obtained for the JR distributions, with 
differences of up to 9%. Comparing to the first experimental measurements by CMS [6] and 
ATLAS [7], all current predictions are compatible within the present experimental errors for all 
channels. 

Due to the reduced scale uncertainty at NNLO and the slightly different NNLO estimates 
of the various groups, cf. Figures 1(a) and (b), and Tables 1 and 2, we conclude that the 
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rates for gauge boson production at LHC energies can be rather confidently predicted with an 
accuracy better than about 9%. These differences are due to the different light sea distributions 
(u, d, s) obtained in the various analyses and require more detailed theoretical and experimental 
investigations in the future. 

3 Higgs Boson Production. The inclusive partonic production cross section 4 of the Higgs 
boson of the Standard Model a{pp — > H° + X) [2] depends on the strong coupling constant oc a 2 s 
and a large part of the cross section exhibits a quadratic dependence of the gluon density. As 
has been shown in Table 1, the a s (M|) values determined, resp. assumed [12], are still different. 
Despite having reached an impressive accuracy of 0(1%) in individual analyses of the world 
deep-inelastic data and related collider data, [8-10,29,31] 5 , and for a large number of other high 
energy processes and decay widths, cf. [33], an overall agreement has not yet been reached. As 
well-known, the gluon density and the value of a s (M§) are anticorrelated, see e.g. the tables of 



M H (GeV) 




\BM10 


[8] 


ABKM09 


[9] 




JR 


[ 


10] 


MSTW08 


[11] 


HERAPDF [12] 


100 


1 


438 


± 





066 


1.380 


± 





076 


1 


593 


± 


0.091 


1.682 


± 





046 


1.417 


110 


1 


051 


± 





052 


1.022 


± 





061 


1 


209 


± 


0.078 


1.265 


± 





038 


1.055 


115 





904 


± 





047 


0.885 


± 





055 


1 


060 


± 


0.072 


1.104 


± 





034 


0.917 


120 





781 


± 





042 


0.770 


± 





050 





933 


± 


0.067 


0.968 


± 





031 


0.800 


125 





677 


± 





038 


0.672 


± 





045 





823 


± 


0.062 


0.851 


± 





029 


0.700 


130 





588 


± 





034 


0.589 


± 





041 





729 


± 


0.058 


0.752 


± 





026 


0.615 


135 





513 


± 





031 


0.518 


± 





037 





647 


± 


0.054 


0.666 


± 





024 


0.541 


140 





449 


± 





028 


0.456 


± 





034 





576 


± 


0.050 


0.591 


± 





022 


0.479 


145 





394 


± 





025 


0.403 


± 





031 





514 


± 


0.047 


0.527 


± 





020 


0.424 


150 





347 


± 





023 


0.358 


± 





028 





461 


± 


0.044 


0.471 


± 





018 


0.377 


155 





306 


± 





020 


0.318 


± 





026 





413 


± 


0.041 


0.421 


± 





017 


0.336 


160 





271 


± 





019 


0.283 


± 





024 





371 


± 


0.039 


0.378 


± 





016 


0.300 


165 





240 


± 





017 


0.253 


± 





022 





335 


± 


0.036 


0.341 


± 





014 


0.269 


170 





213 


± 





015 


0.226 


± 





020 





302 


± 


0.034 


0.307 


± 





013 


0.241 


175 





190 


± 





014 


0.203 


± 





019 





274 


± 


0.032 


0.278 


± 





012 


0.217 


180 





169 


± 





013 


0.182 


± 





017 





248 


± 


0.030 


0.251 


± 





012 


0.195 


185 





151 


± 





012 


0.164 


± 





016 





225 


± 


0.028 


0.228 


± 





011 


0.176 


190 





136 


± 





011 


0.148 


± 





015 





205 


± 


0.027 


0.207 


± 





010 


0.159 


200 





109 


± 





009 


0.121 


± 





013 





170 


± 


0.024 


0.172 


± 





009 


0.131 



Table 3: NNLO predictions for the production cross sections a(pp — > H° + X) \pb\ at \fS = 1.96 TeV. 
The errors refer to the ±la pdf uncertainties. 

the correlation coefficients given in [9]. The values of a s (M§) determined in global analyses 
to derive the parton distribution functions have all lower central values than the world average 

4 The vector boson fusion channel adds a correction of about 10%, with an 0(2%) pdf uncertainty, cf. [30]. 
5 For a compilation see Ref. [32]. 
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Figure 2: Comparison of the NNLO gluon distributions at Q 2 = 4 GeV 2 and Q 2 = (160 GeV) 2 for 
the ratios xg(x, Q 2 )/xg(x, Q 2 )abkm for ABKM09 [9] (full line), JR [10] (dashed line), MSTW08 [11] 
(dotted line), and HERAPDF [12] (dash-dotted line, without error band). 



in [33] 6 , with a spread of Aa s (M§) = 0.0047. If one extended the uncertainty referring to 
the average a s (M|)-value for the r-decay [33] one would obtain Aa s (M|) = 0.0073. Both 
uncertainties imply 

Aa(pp(p) -*H° + X)~ (^^J = 8-4 resp. 13 % . (2) 

Furthermore, there are also still significant differences in the gluon distributions of different 
analyses. The non-perturbative parton distributions at the initial scale Qq of the evolution are 
orthogonal in parameter space to Aq CD , resp. a s (M§), despite correlations. The evolution is 
therefore a direct mapping of the initial conditions and its strength is related to the value of 
a g (M§). The larger this parameter is, the faster the evolution. In Figure 2 we compare the NNLO 
gluon distributions of ABMKM09, JR, HERAPDF, and MSTW08 at a scale Q 2 = 4 GeV 2 and 
Q 2 = (160 GeV) 2 , a typical mass scale for current Higgs boson searches, in the relevant x-range. 
Average values of Bjorken (x) ~ 10 _1 , (x) = y/x\X2 = Mh/VS, in Figure 2(b) correspond to 
the production region at Tevatron, while those at (x) ~ 10~ 2 (3 x 10~ 2 ) are characteristic for 
the LHC at y/S = 14 (7) TeV. At Q 2 = 4 GeV 2 the gluon distributions of JR and MSTW08 
agree rather well in the region x G [10 _2 ,10 _1 ], while ABM10 and HERAPDF yield 20% and 
10% larger values at x = 10~ 2 . All distributions get close for x ~ 10 _1 and for higher values 
the JR distributions becomes largest, followed by MSTW08 and HERAPDF. The evolution to 
fi 2 = (160 GeV) 2 diminishes the differences at x ~ 10 -2 overall to about 5 %. All distributions 
cross around x ~ 3 x 10~ 2 , and at x ~ 10 _1 the gluon distributions by JR (MSTW08) take 

6 For other high energy measurements of a s (M§) yielding lower values see e.g. [34]. 
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~ 20% (10%) larger values than ABKM09, while the HERAPDF values are close to the latter. 
Note that JR obtains a 4 % smaller value of a s {M 2 z ) than MSTW08, and the value of a s (M|) 
determined by ABKM is 3 % smaller than that of MSTW08. 



M H (GeV) 


ABM10 [8] 


ABKM09 [9] 


JR [10] 


MSTW08 [11] 


HERAPDF [12] 


100 


22.82 ± 0.53 


21.18 ± 0.60 


20.48 ± 0.70 


22.95 ± 0.31 


20.90 


110 


18.65 ± 0.44 


17.30 ± 0.49 


16.92 ± 0.56 


18.84 ± 0.26 


17.12 


115 


16.95 ± 0.40 


15.72 ± 0.45 


15.46 ± 0.50 


17.16 ± 0.23 


15.58 


120 


15.45 ± 0.37 


14.34 ± 0.41 


14.17 ± 0.45 


15.69 ± 0.22 


14.22 


125 


14.14 ± 0.35 


13.12 ± 0.38 


13.03 ± 0.41 


14.39 ± 0.20 


13.03 


130 


12.96 ± 0.32 


12.03 ± 0.35 


12.01 ± 0.37 


13.23 ± 0.19 


11.97 


135 


11.92 ± 0.29 


11.07 ± 0.33 


11.10 ± 0.34 


12.20 ± 0.17 


11.02 


140 


10.99 ± 0.27 


10.21 ± 0.31 


10.29 ± 0.32 


11.28 ± 0.16 


10.18 


145 


10.15 ± 0.26 


9.44 


± 


0.29 


9.55 ± 


0.29 


10.45 ± 0.15 


9.42 


150 


9.40 


± 


0.24 


8.75 


± 


0.27 


8.89 ± 


0.27 


9.71 


± 


0.14 


8.74 


155 


8.73 


± 


0.23 


8.13 


± 


0.25 


8.30 ± 


0.25 


9.04 


± 


0.14 


8.13 


160 


8.12 


± 


0.21 


7.56 


± 


0.24 


7.75 ± 


0.24 


8.43 


± 


0.13 


7.57 


165 


7.56 


± 


0.20 


7.05 


± 


0.23 


7.26 ± 


0.23 


7.88 


± 


0.12 


7.07 


170 


7.06 


± 


0.19 


6.59 


± 


0.21 


6.82 ± 


0.21 


7.38 


± 


0.12 


6.62 


175 


6.60 


± 


0.18 


6.17 


± 


0.20 


6.41 ± 


0.20 


6.92 


± 


0.11 


6.20 


180 


6.19 


± 


0.17 


5.79 


± 


0.19 


6.04 ± 


0.19 


6.51 


± 


0.11 


5.83 


185 


5.80 


± 


0.16 


5.43 


± 


0.18 


5.70 ± 


0.18 


6.13 


± 


0.10 


5.48 


190 


5.46 


± 


0.15 


5.11 


± 


0.17 


5.39 ± 


0.18 


5.78 


± 


0.10 


5.16 


200 


4.84 


± 


0.14 


4.55 


± 


0.16 


4.83 ± 


0.16 


5.16 


± 


0.09 


4.60 


220 


3.88 


± 


0.12 


3.67 


± 


0.14 


3.96 ± 


0.14 


4.20 


± 


0.08 


3.73 


240 


3.18 


± 


0.10 


3.02 


± 


0.12 


3.32 ± 


0.13 


3.49 


± 


0.07 


3.09 


260 


2.66 


± 


0.09 


2.55 


± 


0.10 


2.84 ± 


0.12 


2.96 


± 


0.06 


2.61 


280 


2.28 


± 


0.08 


2.19 


± 


0.09 


2.48 ± 


0.11 


2.58 


± 


0.06 


2.26 


300 


2.00 


± 


0.08 


1.94 


± 


0.09 


2.23 ± 


0.11 


2.29 


± 


0.06 


2.00 



Table 4: NNLO predictions for the production cross sections a(pp — > H° + X) \pb] at LHC for = 
7 TeV. The errors refer to the ±lcr pdf uncertainties. 

In Table 3 the predictions for the NNLO Higgs boson production cross section for pp an- 
nihilation at = 1.96 TeV are compared for the distributions [8-12] 7 for the mass range 
100 < M H < 200 GeV. Here the largest cross section differences are those between MSTW08 
and ABKM09 which amount to +22% at M H = 100 GeV and the MSTW08 prediction is +39% 
higher than the one of ABM10 in the present exclusion region around Mh = 160 GeV. This 
difference corresponds to 5.6<x in the pdf-error and is due to both the different gluon densities 
and a s values having been obtained in both analyses, as discussed above. The impact of the 

7 See also [35]. 
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combined HI and ZEUS data is seen by comparing the cross sections for ABM10 and ABKM09 : 
Larger cross sections are obtained at low masses M H < 140 GeV, while for higher masses the 
cross sections become smaller than the ones by ABKM09. The cross sections predicted by JR 
and MSTW08 are compatible within la and those of HERAPDF are about 1.5a larger than 
those of ABM 10. From the above discussion of the gluon distribution and the different values of 
ol s {M\ ) involved, these quantitative relations can be understood to a large extent. To establish 
firm exclusion bounds, cf. [36], the observed variation in the predicted cross sections has to be 
taken into account as an uncertainty. 

In Tables 4 and 5 we compare the different NNLO predictions for pp-scattering at the current 
LHC energy of 7 TeV and the anticipated one of 8 TeV for the next running period [37] for the 
mass range of 100 < Mh < 300 GeV. A quite similar pattern is obtained for both cms energies. 



M H (GeV) 



100 
110 
115 
120 
125 
130 
135 
140 
145 
150 
155 
160 
165 
170 
175 
180 
185 
190 
200 
220 
240 
260 
280 
300 



ABM10 



28.81 ± 0.65 
23.71 ± 0.54 
21.62 ± 0.49 
19.78 ± 0.46 
18.15 ± 0.42 
16.70 ± 0.39 

15.41 ± 0.36 
14.25 ± 0.34 
13.21 ± 0.32 
12.27 ± 0.30 

11.42 ± 0.28 
10.66 ± 0.26 
9.96 ± 0.25 
9.33 ± 0.23 
8.75 ± 0.22 
8.22 ± 0.21 
7.73 ± 0.20 
7.29 ± 0.19 
6.51 ± 0.18 
5.28 ± 0.15 
4.37 ± 0.13 
3.70 ± 0.12 
3.20 ± 0.10 
2.83 ± 0.10 



ABKM09 [9] 



26.81 ± 0.74 
22.04 ± 0.61 
20.09 ± 0.56 
18.38 ± 0.51 
16.86 ± 0.48 
15.52 ± 0.44 
14.32 ± 0.40 
13.24 ± 0.38 
12.28 ± 0.36 
11.41 ± 0.33 
10.63 ± 0.31 
9.92 ± 0.29 
9.28 ± 0.27 
8.69 ± 0.26 
8.15 ± 0.25 
7.67 ± 0.24 
7.22 ± 0.23 
6.81 ± 0.21 
6.09 ± 0.20 
4.96 ± 0.17 
4.13 ± 0.15 
3.51 ± 0.13 
3.05 ± 0.12 
2.72 ± 0.11 



JR [10] 



25.66 ± 0.91 
21.31 ± 0.72 
19.53 ± 0.65 
17.95 ± 0.59 
16.55 ± 0.53 
15.29 ± 0.49 
14.17 ± 0.44 
13.16 ± 0.41 
12.26 ± 0.37 
11.44 ± 0.35 
10.69 ± 0.32 
10.02 ± 0.30 
9.41 ± 0.28 
8.85 ± 0.27 
8.34 ± 0.25 
7.88 ± 0.24 
7.45 ± 0.23 
7.06 ± 0.22 
6.36 ± 0.20 
5.26 ± 0.17 
4.44 ± 0.15 
3.83 ± 0.14 
3.38 ± 0.14 
3.05 ± 0.13 



MSTW08 [11] 



28.85 ± 0.38 
23.83 ± 0.32 
21.77 ± 0.29 
19.96 ± 0.27 
18.35 ± 0.25 
16.93 ± 0.23 
15.65 ± 0.21 
14.51 ± 0.20 
13.48 ± 0.19 
12.55 ± 0.18 
11.71 ± 0.17 
10.96 ± 0.16 
10.27 ± 0.15 
9.64 ± 0.14 
9.06 ± 0.14 
8.54 ± 0.13 
8.06 ± 0.12 
7.62 ± 0.12 
6.84 ± 0.11 
5.61 ± 0.10 
4.70 ± 0.09 
4.03 ± 0.08 
3.53 ± 0.07 
3.17 ± 0.07 



HERAPDF [12] 



26.38 

21.74 

19.85 

18.18 

16.70 

15.39 

14.21 

13.16 

12.22 

11.37 

10.60 

9.90 

9.27 

8.69 

8.17 

7.69 

7.25 

6.85 

6.14 

5.02 

4.19 

3.58 

3.13 

2.79 



Table 5: NNLO predictions for the production cross sections a(pp — > H° + X) \pb] at LHC for y/S = 
8 TeV. The errors refer to the ±la pdf uncertainties. 

At ^/S = 7 TeV the predictions of HERAPDF are lower than those of ABM10, while JR yields 
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also smaller values for masses Mh < 200 GeV and larger ones for higher Higgs masses. MSTW08 
predicts higher cross sections than ABM10, with a tendency of a growing difference towards high 
masses. Overall the predictions are at variance of up to 3a of the pdf-errors, which corresponds 
to maximal deviations of 11-14 %. In the relevant region the gluon densities agree better than 
5 %, but there is still also the uncertainty in a s (Mf ), see Eq. (2). At y/S = 8 TeV the differences 
amount to less than 3.5a of the pdf-error, which corresponds to deviations of 11-16 %. In Table 6 
we compare the NNLO predictions for = 14 TeV for Higgs masses between 100 and 300 GeV. 
They differ between 10-14% and agree better for larger masses. These differences form theory 
errors, which have to be accounted for within feasibility studies, cf. [38], and in searching for 
Higgs boson production at the LHC. 





ABM10 [8] 


ABKM09 


[9J 


JR 


[10] 


MSTW08 




HERAPDF [12] 


1 aa 

100 


71 


16 


± 


1 


53 


67 


27 


± 


1 


78 


62 


24 


± 


2 


62 


70 


73 


± 





AO 

98 


65 


54 


1 in 

110 


60 


05 


± 


1 


27 


56 


60 


± 


1 


A O 

48 


52 


77 


± 


2 


11 


59 


73 


± 





81 


55 


28 


115 


55 


42 


± 


1 


17 


52 


17 


± 


1 


36 


48 


82 


± 


1 


92 


55 


16 


± 





73 


51 


01 


120 


51 


32 


± 


1 


10 


48 


25 


± 


1 


24 


45 


32 


± 


1 


74 


51 


10 


± 





69 


47 


23 


125 


47 


63 


± 


1 


00 


44 


73 


± 


1 


16 


42 


16 


± 


1 


59 


47 


46 


± 





62 


43 


83 


130 


44 


33 


± 





94 


41 


59 


± 


1 


08 


39 


32 


± 


1 


45 


44 


19 


± 





57 


40 


80 


135 


41 


36 


± 





87 


38 


77 


± 


1 


00 


36 


77 


± 


1 


33 


41 


26 


± 





53 


38 


07 


140 


38 


67 


± 





81 


36 


22 


± 





93 


34 


45 


± 


1 


23 


38 


60 


± 





49 


35 


60 


145 


36 


23 


± 





77 


33 


92 


± 





87 


32 


36 


± 


1 


13 


36 


21 


± 





46 


33 


37 


150 


34 


02 


± 





71 


31 


83 


± 





81 


30 


46 


± 


1 


04 


34 


03 


± 





43 


31 


34 


155 


32 


00 


± 





67 


29 


93 


± 





11 


28 


72 


± 





97 


32 


04 


± 





A A 

40 


29 


49 


160 


30 


16 


± 





64 


28 


20 


± 





72 


27 


14 


± 





90 


30 


22 


± 





38 


27 


81 


165 


28 


48 


± 





62 


26 


62 


± 





bo 


25 


70 


± 





83 


28 


58 


± 





3d 


26 


28 


170 


26 


93 


± 





57 


25 


16 


± 





65 


24 


37 


± 





78 


27 


05 


± 





34 


24 


87 


175 


25 


52 


± 





54 


23 


83 


± 





61 


23 


15 


± 





73 


25 


65 


± 





32 


23 


58 


180 


24 


21 


± 





52 


22 


61 


± 





58 


22 


03 


± 





69 


24 


37 


± 





31 


22 


39 


185 


23 


00 


± 





49 


21 


48 


± 





56 


20 


99 


± 





64 


23 


18 


± 





29 


21 


30 


190 


21 


90 


± 





47 


20 


44 


± 





53 


20 


04 


± 





61 


22 


09 


± 





28 


20 


29 


200 


19 


91 


± 





43 


18 


59 


± 





49 


18 


33 


± 





55 


20 


14 


± 





26 


18 


49 


220 


16 


75 


± 





37 


15 


64 


± 





41 


15 


59 


± 





45 


17 


03 


± 





22 


15 


62 


240 


14 


38 


± 





33 


13 


44 


± 





36 


13 


54 


± 





38 


14 


70 


± 





19 


13 


46 


260 


12 


60 


± 





29 


11 


79 


± 





33 


12 


01 


± 





34 


12 


96 


± 





18 


11 


85 


280 


11 


27 


± 





27 


10 


56 


± 





30 


10 


86 


± 





30 


11 


66 


± 





17 


10 


64 


300 


10 


32 


± 





25 


9. 


59 




0. 


28 


10 


07 


± 





29 


10 


75 


± 





16 


9. 


30 



Table 6: NNLO predictions for the production cross sections a(pp — > H° + X) \pb] at LHC for = 
14 TeV. The errors refer to the ±1<7 pdf uncertainties. 

In Figure 3 we illustrate the different predictions of ABM10, JR, HERAPDF, and MSTW08 
for the inclusive Higgs production cross section at the LHC for the energies y/S = 7, 10 and 
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14 TeV. The cross sections rise roughly oc y~S. In this energy range ABM10 and MSTW08 
predict nearly equal Higgs boson production cross sections, while those by JR and HERAPDF 
are found to be lower. For ABM 10 and JR we include in the error also the scale variation 
uncertainty varying = //p in the range [Mh/2, 2Mh] and added the pdf-errors in quadrature. 
It should be kept in mind that the scale variation errors are correlated when comparing the 
predictions for the different pdf sets. Although the various NNLO predictions differ w.r.t. the 
present pdf-errors, they are consistent within la taking the scale variation errors into account. At 
NNLO these are still significant, which will require to account for even higher order corrections, 
despite the well-known fact that soft gluon resummation leads to improvements (see [2] and 
references therein). 



40 



30 



10 



NNLO 

ABM10 

JR 

HERAPDF 
MSTW08 



M H 
120 GeV 



150 GeV |± A i 
180 GeV 



T 

I' 

V 



7 TeV 



10 TeV 



14 TeV 



Figure 3: Predictions of the inclusive Higgs-boson production cross sections at NNLO for different 
energies at the LHC for the parton distributions ABM10, JR, HERAPDF, MSTW08, [8,10-12]. For 
the ABM10 and JR distributions the scale variation errors corresponding to the range Mjj/2 < fj,p = 
/xr < 2Mh are included. The inner error bars refer to the pdf-errors only. 

The typical scale uncertainty of each individual NNLO prediction amounts to about ±9% at 
7 TeV and ±8% at 14 TeV which doubles at NLO [10,39]. A measure for the uncertainty of the 
Higgs boson production cross sections is thus obtained adding in quadrature the largest difference 
of central values and the largest error, including scale variations, of one of the predictions, and 
dividing this result by the central value of the smallest prediction. Higgs boson production at 
LHC can be predicted with an accuracy of about 10-17 % at NNLO 8 (with a total uncertainty 
being almost twice as large at NLO), whereas the uncertainty almost doubles at the Tevatron 
(y/S = 1.96 TeV). Furthermore, the NNLO predictions are typically about 20 % larger than at 
NLO. 

We finally remark that a possible way to determine experimentally the gluon density in the 
mass region relevant for Higgs boson production consists in a precise measurement of the single 

8 Notice that the largest uncertainty refers to the predictions for Mh = 120 GeV at 14 TeV, and the uncertainty 
decreases for Mh increasing; the smallest uncertainty of about 10 % refers to the predictions for Mh = 180 GeV 
at 7 TeV. 
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top quark production cross section [40-42] for scales around /j, ~ m t . At present somewhat 
differing results for this process were obtained at Tevatron [43], however, still with large errors, 



a(pp -> t{t) + b(b) + X) = 2.3 



+0.6 
- 0.5 



(stat. + sys.) pb CDF 



(3) 



= 4.19 



+1.24 
- 1.14 



(stat. + sys.) pb DO 



(4) 



with 9 V tb = 0.91 ±0.11. Theoretical NLO QCD predictions 10 are 



a(pp ->• t(t) + b{b) + X) 



14| 2 (2.420 ± 0.003) pb MSTW08 
V t6 | 2 (2.195 ± 0.003) pb ABKM09 , 



(5) 
(6) 



which differ by about 10%. The NLO scale variation errors according to m t /2 < /j-f,r < 2mt 
amount up to 10% at Tevatron, are lower than 5 % at LHC [42], and can be further improved at 
NNLO. The different values in (5) and (6) are caused partly by yet different gluon distributions 
in the region x ~ 0.1, see Figure 2, but at least to the same extent they are due to the different 
values of a s , cf. Table 1. Precision measurements of the single top production cross sections at 
the LHC, although challenging, with an accuracy of 5% or better would allow to determine the 
gluon density in the same experiments which measure the production of the Higgs boson as well. 

4 Conclusions. The present analysis shows that the current NNLO analyses of the world 
deep- inelastic data, supplemented with other relevant hard scattering data to measure the parton 
distributions, yield still different results at the level of less than about 10 % for and Z° 
boson production, and at the level of 10 — 17% for Higgs boson production at the LHC. (These 
uncertainties would almost double at NLO.) The present variations in the predictions based 
on [8-13] are both due to the QCD-scale Aq C d (or a s (M§)) and differences in the parton densties, 
as well as due to the renormalization and factorization scale uncertainties of the theoretical 
predictions. Future detailed work is needed to study and to further delineate these differences. 
The differences documented in this study form essential contributions to the theory error of the 
measurements of these scattering cross sections at the Tevatron and the LHC, which should 
therefore be compared to all the predictions based on the distributions ABM10 [8], JR [10], 
MSTW08 [11], and HERAPDF [12] at NNLO. 
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Appendix. In the following we summarize the different experimental measurements of the 
W ± and Z°-boson production cross sections in pp and pp collisions for comparison with the 
predictions based on the current NNLO parton densities. 

The cross sections for pp collisions measured by UA1, UA2 and DO experiments are (averaging 

9 Assuming unitarity for the CKM matrix would imply a far smaller error [25]. 

10 We would like to thank S. Alioli for computing the reference values for single top production with POWEG [41]. 
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over the different channels measured by UAl for Z° production) 



UA1 (0.546 TeV): 


a w++w- B ( W ± 


= 0.55 ±0.12 nb 


(A.l) 




a z °B(Z° -+ ) 


= 0.070 ±0.046 nb 


(A.2) 


UA2(0.546 TeV): 




= 0.57 ±0.08 nb 


(A.3) 




a z° B (Z° -+ t + n 


= 0.116 ±0.041 nb 


(A.4) 


UAl (0.630 TeV): 


a w++w ~B(W ± ^£u) 


= 0.63 ±0.13 nb 


(A.5) 




a z °B(Z° -> £ + r) 


= 0.070 ±0.019 nb 


(A.6) 


UA2(0.630 TeV): 


a w++w- B ( W ± _^ lv) 


= 0.57 ±0.08 nb 


(A.7) 




a z °B(z° -> e+r) 


= 0.073 ±0.016 nb 


(A.8) 


D0(0.630 TeV): 


a w++w ~B(W ± ^£u) 


= 0.658 ± 0.058 ± 0.034 nb . 


(A.9) 



The cross section measurement of CDF at v^S* = 0.630 TeV has not been published [44]. 

At similar energies of v^S* = 0.5 TeV the inclusive W ± boson cross sections were recently 
measured by the PHENIX experiment [17] in pp collisions at RHIC : 

PHENIX(0.5 TeV): 



a w+ B(W + -)• e + u) = 0.1441 ± 0.0212 0103 ± 15% ( norm ) nb ( A - 10 ) 
<r w ~ B(W~ -> e~V) = 0.0317 ±0.0212 ^qJ^ ± 15%(norm) nb . (A.ll) 

Here and in the following the errors are given in the sequence of statistical, systematic, and 
luminosity/normalization errors. 

The pp experiments CDF and DO at the Tevatron performed inclusive measurements for W ± 
and Z Q production at \fS = 1.8 and 1.96 TeV : 

CDF(1.8 TeV) [18] : 

a w++w ~B(W ± ->■ Iv) = 2.49 ±0.12 nb (A.12) 

a z °B(Z° ->■ £ + r) = 0.231 ±0.012 nb (A.13) 

DO (1.8 TeV) [20] : 

a w++w ~B(W ± ->• iv) = 2.310 ±0.001 ±0.005 ±0.100 nb (A.14) 

a z °B(Z° ->■ = 0.221 ±0.003 ±0.004 ±0.010 nb (A.15) 



CDF(1.96 TeV) [18,19] : 



a w++w ~B(W ± -*lv) = 2.749 ±0.010 ±0.053 ±0.165 nb (A.16) 
a z °B(Z° ->■ £ + r ) = 0.2549 ±0.0033 ±0.0045 ±0.0152 nb (A.17) 
= 0.2566 ±0.0007 ±0.0020 ±0.0154 nb. (A.18) 
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First measurements of the inclusive W ± and Z° boson cross sections at LHC were performed 
at t/S = 7 TeV : 

CMS [6] : a w++w ~B(W ± ->■ tv) = 9.951 ± 0.073 ± 0.280 ± 1.095 nb (A.19) 

a w+ B(W + ->■ tv) = 5.859 ±0.059 ±0.168 ±0.645 nb (A.20) 

a w ~B(W~ (TV) = 4.140 ±0.064 ±0.254 ±0.455 nb (A.21) 

a z °B(Z° ^ t + £~) = 0.960 ±0.037 ±0.059 ±0.106 nb (A.22) 

ATLAS [7] : a w++w ~B(W ± -> &/) = 9.96 ± 0.23 ± 0.50 ± 1.10 nb (A.23) 

<7 W+ B(W + ->■ tv) = 5.93 ±0.17 ±0.30 ±0.65 nb (A.24) 

a w ~ B(W~ -> TV) = 4.00 ±0.15 ±0.20 ±0.44 nb (A.25) 

a 20 B(Z° ->£+£-) = 0.82 ±0.06 ±0.05 ±0.09 nb. (A.26) 

These cross sections correspond to analyzed samples of 2.88 pb~ x and 0.32 pb~ x , while the current 
run may accumulate O(50 pb' 1 ), which will lead to a significant reduction of the statistical error 
and allow to improve the systematic errors. 

To compare with the theoretical predictions one may use the current measured branching 
fractions [25] 

B(W + ->■ e + v e ) = 0.1075 ±0.0013 , B{W + -> n + v„) = 0.1057 ± 0.0015 , 

B(W + ->■ r+z/ T ) = 0.1125 ± 0.0020 , (A.27) 

with ->■ £ _ I70 = ->■ and 

B(Z°^e + e-) = 0.03363 ± 0.00004 , B(Z° ->■ ^V") = 0.03366 ± 0.00007 , 
5(Z°^t + t-) = 0.03367 ± 0.00008 , (A.28) 

or their corresponding averages. 
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